Decoherence and the Loschmidt echo 
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Environment-induced decoherence causes entropy increase. It can be quantified using, e.g., the 
purity c = Trp 2 . When the Hamiltonian of a quantum system is perturbed, its sensitivity to such 
perturbation can be measured by the Loschmidt echo M(t). It is given by the average squared overlap 
between the perturbed and unperturbed state. We describe the relation between the temporal 
behavior of c(t) and M(t). In this way we show that the decay of the Loschmidt echo can be 
analyzed using tools developed in the study of decoherence. In particular, for systems with a 
classically chaotic Hamiltonian the decay of ? and M has a regime where it is dominated by the 
classical Lyapunov exponents. 



Environment-induced decoherence is an essential in- 
gredient of the quantum-classical transition [1| . Its im- 
plications for the quantum versions of classically chaotic 
systems are especially intriguing, as they go beyond 
the restoration of the quantum-classical correspondence. 
Two of us discussed this issue for the first time in 0, 
presenting a surprising result that has been since amply 
corroborated 0,0, OJ H : For a quantum system with a 
classically chaotic Hamiltonian the rate at which infor- 
mation about the initial state is degraded by the environ- 
ment becomes independent of the system-environment 
coupling strength. This rate (e.g., the von Neumann en- 
tropy production rate computed from the reduced density 
matrix of the system) is set by the classical Lyapunov ex- 
ponents, provided that the coupling strength are within 
a certain (wide) range. This result has important impli- 
cations and can be a way to define quantum chaos Q ■ A 
related but independent way to do this was considered 
first by Asher Peres by appealing to a property of 
quantum systems with a classically chaotic counterpart: 
Although quantum dynamics is insensitive to small dif- 
ferences in initial conditions, it seems to be highly sen- 
sitive to perturbations in the Hamiltonian. The similar 
idea of "hypersensitivity" was studied in an information 
theoretic framework by Caves and coworkers 0- More 
recently, Levstein, Pastawski and collaborators [lfl III] ] 
experimentally studied sensitivity to perturbations by 
measuring the Loschmidt echo in a many-body spin sys- 
tem. Their work underscored the need for a systematic 
study of various aspects of the true nature of the sensitiv- 
ity displayed by quantum systems when their Hamilto- 
nian is perturbed. Further analytical |l^ and numerical 
[l^j ] work triggered an intense activity that helped clar- 
ify some aspects of the temporal dependence of the echo 
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The measure of the echo signal is the overlap between 
two states that evolve from the same initial wave function 
under the influence of two Hamiltonians (the unper- 
turbed one, H , and the perturbed one Ha = H + A). 
More precisely, when Uq and Ua denote the correspond- 



ing evolution operators, the echo is defined as 
M A (i) - |(* |[/ A (i)E/ (*)|*o)| 2 . 



(1) 



The quantity M(t), obtained by averaging Ma over an 
ensemble of perturbations, can be studied analytically 
and displays a rich temporal dependence. One interest- 
ing regime was analyzed by Jalabert and Pastawski |l^ 
who showed, using a semiclassical approximation, that 
there is a window of values for the perturbation strength 
for which M(t) decays with a rate equal to the classical 
Lyapunov exponent. In spite of the simple discussion pre- 
sented above, the physically relevant evolution will typ- 
ically not be unitary: Decoherence caused by the inter- 
action with the environment will suppress the echo even 
in the absence of the perturbation A. We shall, however, 
adhere to the usual assumption 

that the evolutions are unitary, and show that even in 
that case of decoherence-free echo supression it is possi- 
ble to draw useful conclusions from the analogy we dis- 
cuss below. 

In this letter we establish a direct connection between 
decoherence and the decay of the Loschmidt echo. In 
particular, we relate the evolution of M(t) and the linear 
entropy (or purity) of an open quantum system. The ex- 
istence of a kinship between decoherence and the decay 
of the Loschmidt echo is not unexpected: such possibility 
was noted, for example, in Refs. [ill Il2j. but was never 
formally established. Demonstrating the relation of these 
two important areas is interesting not just from a funda- 
mental point of view but may also prove useful since it 
allows one to use results obtained in the theory of open 
quantum systems to understand better the behavior of 
the echo. 

The key step in our demonstration is a simple obser- 
vation: The average echo M(t), for an ensemble of per- 
turbations characterized by a probability density -P(A), 
is 

M(t) =JDA P(A) |<*o|C/i(i)[/o(i)|*o)| 2 . (2) 

This equation can be rewritten in a more convenient way 
by defining the density matrix of the average perturbed 
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state: 



p(t) = / DA P(A) U- A («)|*o)(*o|tl(t). (3) 



In fact, M(t) is simply the overlap between the av- 
erage state p(t) and the unperturbed density ma- 
trix po(t) evolved from the initial state, po{t) = 
U (t)\V )(y \tf (t): 



M(t) =Tr (p(t)p (t)) 



(4) 



Once we recognize this simple fact we can go one step fur- 
ther and notice that the state p generally evolves in time 
according to a master equation which is of the same kind 
as the ones arising in the study of decoherence. Hence, 
the evolution of the echo M(t) is directly placed in the 
context of open quantum systems and decoherence. 

Equation (@J can be used to establish an inequal- 
ity between the average echo M(t) and the purity 
?(t) = Tip 2 (t) (used to characterize decoherence). Us- 
ing Schwartz inequality (applied to the inner product 
between two hermitian operators) and assuming that the 
initial state is pure we find that 



M 2 (t) < 



(5) 



Related inequalities were noticed and used in a some- 
what different context in 0, 0] . Eq. (0) implies that 
when the purity Tip 2 decays exponentially with a rate 
7_d, then the overlap M(t) should also decay exponen- 
tially (or faster) with a rate that should be no less than 
7_d/2. However, as we will see later (and as has been 
established in the literature 0,0,E1)j there is an im- 
portant regime (the so-called Lyapunov regime) where 
both quantities decay with the same rate set by the Lya- 
punov exponent. 

Let us now analyze some generic features of the evolu- 
tion of the average state p(t). It is useful to notice that 
p generally obeys a master equation with non-unitary 
terms. These terms arise because averaging of the evo- 
lution over an ensemble of perturbations yields an effect 
analogous (although not equivalent) to the tracing out 
of the unobserved degrees of freedom pj. We will find 
it convenient to consider a simple form of perturbation 
(even though results are not strongly dependent on it). 
Let us assume that A(x,t) = V{x)J{t), where V(x) is a 
function of the coordinates of our system and J(t) is an 
external source. For this case, averaging over the pertur- 
bation A consists of averaging over functions J(t). We 
will also assume that the probability density P(J) is a 
Gaussian functional whose width defines the temporal 
correlation function for the sources. Thus, we will use 



P(J) = Nexp 



dtdt'Jifiv- 1 ^') J(t'] 



where N is a normalization factor. The correlation func- 
tion of the source is then J DJ P{J)J(t)J(t') = v(t,t'). 
Using this, we can show that the evolution operator for 
p(t) has a pat h integral representation with an influence 
functional 18] given by 



F[x,x'] = exp(-- 



dtdt'V-(t)u{t,t')V-(t')), (7) 



where V-(t) = V(x(t)) — V(x'(t)). In some simple but 
physically relevant cases it is possible to write a master 
equation for p(t). In fact, when the noise is white, i.e. 
u(t, t') = 2D5(t — t'), we can show that 



P = -7T [H ,p] 
in 



D[V(x),[V(x),p}}. 



(8) 



While the first term on the rhs of this equation gener- 
ates unitary evolution, the second term is responsible for 
decoherence: It induces a tendency towards diagonaliza- 
tion in position basis and, in the Wigner representation, 
it gives rise to a diffusion term. For the simplest case of 
V(x) = x the equation for the Wigner function reads 



W(x,p) = {H ,W} 



MB 



D d 2 pp W(x,p), 



(9) 



(6) 



where the bracket in the rhs is the so-called Moyal 
bracket, responsible for unitary evolution 0. 

Equations like (JSJ and © arise if we consider a quan- 
tum system interacting with a quantum environment 
formed by a collection of harmonic oscillators [l^. I n 
such a case the absolute value of the influence functional 
generated by the environment is identical to Q provided 
one chooses the spectral density and the initial state of 
the environment in such a way that the noise-kernel it 
produces is equal to the kernel v(t,t') appearing in (7J|. 
However, in general the influence functional is a complex 
number whose phase is responsible for dissipation (noise 
and dissipation kernels are connected as mandated by 
the fluctuation-dissipation theorem). There is a physi- 
cally relevant limit (usually associated with high temper- 
atures) where relaxation effects can be ignored. This is 
the most interesting limit in decoherence studies aimed at 
understanding the quantum-classical correspondence Q] . 
Thus, in such a limit, the evolution of the average state 
p is identical to that of a quantum system that interacts 
with an environment. 

A convenient way to visualize the transition from 
quantum to classical is provided by the Wigner func- 
tion, whose oscillations are the signature of quantum 
interference. They should be suppressed by the de- 
coherence term to make the quantum-classical corre- 
spondence possible. This is indeed what happens: 
When the Wigner function oscillates with a well de- 
fined wave-vector k p along the momentum direction 
(W(x,p,t) ~ A(x,t) cos {kppj), the decoherence term 
in ijnj washes out oscillations exponentially fast with a 
rate Yd — Dk 2 . We can see the behavior of a typical 
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FIG. 1: Wigner function of an initially Gaussian state evolved 
with a chaotic Hamiltonian without (top) and with (bottom) 
decoherence. The system is a particle moving in a driven dou- 
ble well potential (see 5]). A region of area Ti is shown in the 
top panel where sub-Plan ckian structure is evident |2(i|. The 
color scale is positive from yellow to red, shades of blue are 
negative and white is zero. In the top panel we can appreciate 
the distinct regions Ao {Ac) where the Wigner function Wo 
oscillates rapidly (is positive), used in Eq. II 2> . 

Wigner function for a chaotic system (a driven double 
well analyzed in 0]) with and without decoherence in 
Figure 1. Taking into account our previous discussion, 
the echo M[t) is obtained by computing the overlap be- 
tween the two Wigner functions displayed in the figures. 
The purity c should be computed by taking the overlap 
of the decohered Wigner function with itself. Below, we 
will discuss the relation between these two quantities. 

The master equation 10 can be used to obtain the time 
derivatives of the purity q and the echo M : 

4 = 2D J dxd P W{x,p)d 2 p W(x lP ), (10) 

M = D J dxdpW {x,p)d^W{x,p). (11) 

Equation (|1C)|> has been used before to show the exis- 
tence of a domain of exponential decay for the purity Trp 2 
0,0,0,0] with a rate that, for classically chaotic systems, 



is independent of the diffusion D. The central piece of 
the argument is the following: After integrating by parts, 
equation l|10(l can be rewritten as <;/<; = —2D /a 2 , where 
a characterizes the dominant wavelength in the spectrum 
of the Wigner function (i.e. a~ 2 = f(d p W) 2 / JW 2 ). 
Thus, the rate of change of the purity becomes indepen- 
dent of the diffusion constant when a 2 is proportional to 
D. This happens indeed as a consequence of the compe- 
tition between two effects. The first one is the tendency 
of chaotic evolution to generate (exponentially fast, at a 
rate set by the Lyapunov exponent A) small scale struc- 
ture in the Wigner function. The second effect is due 
to diffusion, which tends to wash out small scales expo- 
nentially fast at a rate determined by the product Dk 2 . 
As demonstrated in 0] , these two effects reach a balance 
when the scale is such that a 2 = 2D/X. When this equi- 
librium is reached, the purity <; decreases exponentially 
at a rate fixed by A. 

For this behavior to take place the diffusion constant 
should be above a threshold Q. Otherwise the critical 
width is not established (indeed, in the previous argu- 
ment the implicit assumption is that the time scale for 
diffusion to wash out a fc p -oscillation is shorter than the 
time scale for the oscillations to be regenerated by the 
dynamics). This simple scenario enables us to under- 
stand why there is a regime where purity decreases ex- 
ponentially with a Lyapunov rate. The same argument 
can also be used to analyze the exponential decay of the 
Loschmidt echo. In fact, equations (|1(J|) and (|llfl just 
differ by a factor of 2 and by the presence of Wo instead 
of W inside the integral. As before, we can transform 
the evolution equation of the echo into M/M — —D/a 2 , 
(7 2 = J W d pp W/ JW W. When decoherence is effec- 
tive and the dominant structure in W approaches the 
critical value, the smallest scales of the pure Wigner func- 
tion Wo continue contracting and developing smaller and 
smaller scales (sub-Planck scales are reached quickly in 
chaotic quantum systems [20J). In such a case, one can 
easily conclude that a 2 = 2a 2 . This is readily seen even 
in the crude approximation of W ~ exp(— p 2 /2a 2 ) and 
Wo ~ exp(— p 2 /2a 2 )i with a ~ exp(— Xt) and t > 1/A. 
Therefore, when the purity starts decaying at the Lya- 
punov rate the echo does precisely the same. 

Using the above ideas we are now in a position to 
present a more complete and illustrative picture of the 
time dependence of the echo M(t) and the purity q(t). 
For the sake of simplicity we focus on the echo but the 
same reasoning applies to the purity. To compute the 
overlap M = J dxdpWoW we can split the phase space 
integral into two regions: the region Ac close to the clas- 
sical unstable manifold of the initial state, where Wo is 
positive, and the region Ao over which Wo oscillates (see 
Figure 1): 

M(t)= I dxdpW W+ [ dxdpWoW. (12) 

J Ao J A c 
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In the oscillatory region we can estimate the value of the 
integral assuming that there is a dominant wave vector 
k p . In such a case, from Eq. 10 we can assume that 

W ~ Woe ~ Dk p t . If more than one scale is present the 
result would be a sum of terms like this one. For the sec- 
ond integral, we can also use a crude estimate assuming 
that Wq and W are constant over their respective effec- 
tive support. In particular, W ~ l/Ac since its integral 
over Aq cancels out. As W approaches the critical width 
a along the stable manifold, the area of its effective sup- 
port grows exponentially. Therefore, one gets that the 
second integral is f A W W ~ W A C W ~ W ~ e~ At . 
Thus, combining the two results we find that the ex- 
pected behavior of the Loschmidt echo is 

M(t) = a exp(-Xt) + b exp {-Dk 2 p t) (13) 

for appropriate prefactors a and b. 

This result was previously derived for the Loschmidt 
echo using semiclassical techniques 0. The first term 
gives the Lyapunov decay, while the second one describes 
the so-called Fermi golden rule regime (FGR) [13. In 
this case the rate is proportional to the diffusion coef- 
ficient (which is itself proportional to the square of the 
strength of the perturbation). As mentioned above, a 
similar result is expected for the purity. 

Our treatment is valid in a semiclassical regime where 
the evolution of the Wigner function is dominated by the 
classical Hamiltonian flow and the corresponding inter- 
ference fringes generated when its phase space support 
folds. The virtue of this analysis, entirely based on prop- 
erties of the evolution of W derived in the context of de- 
coherence studies, is not only its simplicity but also the 
fact that it enables us to identify the regions of phase 
space that can be associated with each of the terms ap- 
pearing in (|T31) : the FGR contribut ion arises from the 
decay of the interference fringes while the Lyapunov con- 
tribution is associated with the behavior of W near the 
classical unstable manifold. Recent studies using semi- 
classical techniques point in the same direction |21| . 

It is interesting to perform a better estimate of the 
integral over the region Aq. Assuming that the local 
Lyapunov exponent is constant along the unstable man- 
ifold, one can approximate the value of the integral by 
using the corresponding result for the simplest system 
with an unstable fixed point: the inverted oscillator (10) 
with Hamiltonian Ho = p 2 /2m — mX 2 x 2 /2. In such a 
case, the echo can be computed exactly and turns out to 
be 

M IO (t) = (1 + r sh(2A<) + r 2 (sh 2 {Xt) - X 2 t 2 ))~ 1/2 . 

(14) 

Here r = a 2 /4a 2 , where <jj is the momentum dispersion 
of the initial state. This exact result shows that for long 
times (At ^> log(r)/2) the echo Mjo always decays as 
exp(— Xt). For initial times a decay with a rate deter- 
mined by diffusion is observed but this transitory regime 



always leads to a decay dominated by the Lyapunov ex- 
ponent. The initial transient is sensitive to the details 
of the noise statistics. For example, we can also exactly 
evaluate the echo when the noise kernel is flat (i.e., v[t, t') 
independent of t and t'). For such extreme case the long 
time behavior of the echo is not changed but the initial 
transient displays a quadratic decay. 

We expect the analogy between Loschmidt echo and 
decoherence not only to enable intuitive derivations like 
the one leading us to equation but also to provide 
new insights into theoretically unexplained experimental 
features such as the quadratic decay observed in |ll| . Our 
results are also relevant for quantum computation as the 
Loschmidt echo is a measure of the fidelity with which a 
given algorithm is implemented. The Lyapunov decay of 
the fidelity could hinder the practical implementation of 
such computers, which would then have to deal with an 
exponential increase of error probability at a rate which 
is independent of the coupling to the environment and is 
solely fixed by the (possibly chaotic) nature of the un- 
derlying physical system. The high sensitivity of chaotic 
systems to perturbations seems to be connected to the ef- 
ficiency of these systems to produce decoherence on other 
systems, which has been under recent investigation [22j . 
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